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ABSTRACT 



The MIPSGAL 24 jim Galactic Plane Survey has revealed more than 400 compact-extended 
objects. Less than 15% of these MIPSGAL bubbles (MBs) are known and identified as evolved 
stars. We present Spitzer observations of 4 MBs obtained with the InfraRed Spectrograph to 
determine the origin of the mid-IR emission. We model the mid-IR gas lines and the dust emission 
to infer physical conditions within the MBs and consequently their nature. Two MBs show a 
dust-poor spectrum dominated by highly ionized gas lines of [O IV] , [Ne III] , [Ne V] , [S III] and 
[S IV]. We identify them as planetary nebulae with a density of a few 10'^ cm^"^ and a central 
white dwarf of > 200, 000 K. The mid-IR emission of the two other MBs is dominated by a dust 
continuum and lower-excitation lines. Both of them show a central source in the near-IR (2MASS 
and IRAC) broadband images. The first dust-rich MB matches a Wolf-Rayet star of ^ 60, 000 
K at 7.5 kpc with dust components of ^ 170 and ^ 1750 K. Its mass is about 10"'^ M0 and 
its mass loss is about 10~^ Mo/yr. The second dust-rich MB has recently been suggested as a 
Be/B[e]/LBV candidate. The gas lines of [Fe II] as well as hot continuum components 300 and 
~ 1250 K) arise from the inside of the MB while its outer shell emits a colder dust component 
(~ 75 K). The distance to the MB remains highly uncertain. Its mass is about 10"'^ Mq and its 
mass loss is about 10^^ Mq/jt. 



1. Introduction 

Evolved stars are significant contributors to the 
ionizing photons budget and kinetic energy re- 
leased in the interstellar m edium through their 
intense winds and radiation (jCrowther fc Dessart 
19981 ). They are also the main sources of dust in 



the Galaxy and beyond. In the Milky Way and 
in the local universe, stars from the asymptotic 
giant branch (AGB) are the main producers of 
dust while red supergiants, novae and Wolf-Rayet 
stars (WR) contr i bute only traces of dust(e.g. 
Tielens et ai]|2005HMatsuura et al.ll2009t ). Super- 
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novae (SNe) have also been suggested as signifi- 
cant contributors to the dust production, at least 
at higher redshift, as large ai nount of dust are de- 
tected in young g alaxies (e.g. Bertoldi et al.ll2003 : 
Dwek et al.l 120091) . However, the dust yield from 



most sou rces, SNe in par ticular, is still controver- 
sial (e.g. Rho et al. 20091 and references therein). 



Evolved stars believe d to be present in the Galaxy 
are still "missing " (IShara et all Il999l ). Indeed, 
about 2500 planetary nebulae (PNe) are known in 
our Galaxy but only 224 central stars with spectral 
types are reported in the Strasbourg-E SO Cata- 
logue of Galactic PNe (jAcker et al.lll992l) . At least 
1000 WRs stars are expected to be located within 
our Galaxy, but only 300 have been observed so 



far (|van der Huchtll200ll l200d ). Less than 50 lu- 



minous blue variabl es (LBVs) are kno wn or identi 



fied as candidates (jClark et alj 120051 ) . Therefore, 



it is essential to hunt for these evolved stars in 
order to determine their true nature and obtain 
an accurate census of the dust production in the 
Milky Way as well as in other galactic systems. 

izuno et al.l I (|2010l) discovered over 400 small 
(< 1') rings, bubbles, disks and shells from vi- 
sual ins pection of the M IPSGAL 24 mosaic 
images (jCarev et al.l2009l) . These MIPSGAL bub- 
bles (MBs) are pervasive through the entire Galac- 
tic plane in the mid-infrared. These objects span 
a large range of morphologies, sizes and fluxes. 
Their distribution is approximately uniform in 
Galactic latitude and longitude, and the aver- 
age density is around 1.5 MBs per square de- 
gree. Only 15% of the MBs exist in available cat- 
alogs from the VizieR and SIMBAD database. A 
large majority (8 out of 10) of the already known 
MBs are found in the MASH Catalogue of PNe 
( Parker et al.ll2006|) and i n the Catalogue of Galac- 
tic PNe (jKohoutek 2001 ). Several supernova rem- 
nants (SNRs), WR stars, LBVs and emission line 
stars are also identified from published catalogues. 

Extended emission in rings, shells and discs 
may be found (1) around main sequence and sub- 
giant stars in the form of debris disc/rings as- 
sociated with protoplanetary system, (2) around 
LBVs or associated with SNRs in the form of 
larger ring-type structure and (3) around PNe 
and symbiotic s tars in the form of well dene d 
ring structures ( Phillips fc Ramos-LariosI 2008f) . 
It seems unlikely we can resolve structure around 
main sequence stars in the Galactic plane. The 
usual suspects for the MBs are thus mainly associ- 
ated with stars in the late stages of their evolution. 
In these cases, the extended emission would arise 
from hot, small dust grains and/or from hot ion- 
ized gas in the stellar winds and ejecta. The MBs 
could account for the "missing" massive evolved 
stars. Because 85% of these objects are yet un- 



known, there is a great potential for the discovery 
of several tens of new candidates. The identifica- 
tion of a statistical sample of these shells is thus 
of major importance for the study of stellar evo- 
lution, stellar mass-lo ss, dust life-cycle a n d ISM 



dynamics. R e cently , Gvaramadze et al. ( 2010f ) 



Wachter et al.l (|2010[ ) observed the central source 



detected in some MBs 24 ^m images with optical 
and/or near-IR spectrometer and identify several 
new WR/LBV candidates. However, such obser- 
vations can only be performed on MBs with a cen- 
tral source and does not answer the question on 
the origin of the MIPS 24 /im emission. 

In this paper, we present mid-IR spectroscopic 
observations of 4 MBs obtained with the high res- 
olution module (R ^ 600) of the InfraRed Spec- 
trograph (IRS, iHouck et a l.' '2004') on board the 
Spitzer Space Telescope t Werne r et al.. .2003) as 
part of our Spitzer Director's Discretionary Time 
(DDT) observational program. It is a unique set 
of observationf0 to characterize the origin of the 
MIPS 24 /im emission of the MBs. Indeed, ob- 
servations from the ground are limited in terms 
of sensitivity due to the extended nature of the 
MBs and the atmospheric transmission at those 
wavelengths, while from space, the next satellite 
with mid-IR capabilities (the James Webb Space 
Telescope) is not supposed to be launched before 
2014. Meanwhile, the Stratospheric Observatory 
for Infrared Astronomy (SOFIA) is the best fa- 
cility that can provide mid-IR capabilities, albeit 
with lower sensitivity than Spitzer. The structure 
of the paper is as follows. We present the obser- 
vations and their reduction in section [2] and [3] re- 
spectively. We then analyze each MB separately, 
starting with the dust-poor MBs, in section|31 We 
give our conclusions in section [5] 

2. Observations 

The four MBs we present here have been ob- 
served on the 11th of April 2009, about a month 
before the end of Spitzer cryogenic mission. The 
four MBs are selected among the unidentified ob- 
jects of the catalog that are visible at this time 
of the year by Spitzer. To make the best use of 
the one hour of observations we were granted, we 



^11 additional MBs have been observed in a different pro- 
gram with the low resolution module and will be presented 
in a subsequent paper. 
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Table 1: List of the observed MBs. Nomenclature is from Mizuno et al. ( 201Cll ). 





Name 




RA (J2000) 


DEC(J2000) 


hi (Jy) 


r2i (") 


ho (Jy) 


MB4001 


MGE314.5619^ 


-00.1985 


217.005 


-60.4739 


0.36±0.01 


15 


< 0.2 


MB4006 


MGE316.5508^ 


-00.5816 


220.366 


-59.3494 


0.98±0.02 


21 


< 0.3 


MB3957 


MGE305.6514^ 


-00.3495 


198.737 


-62.3980 


0.59±0.02 


19 




MB4121 


MGE337.55434 


-00.2200 


249.178 


-46.9391 


7.4±0.4 


20 


13 



(a) 



(b) 



(c) 



(d) 

Fig. 1.— Three color image (blue is IRAC 4.5, green is IRAC 8.0 and red is MIPS 24 /xm, 4.2'x4.2') of (a) 
MB4001, (b) MB4006, (c) MB3957, (d) MB4121. 



chose four MBs that exhibit a significant MIPS 
24 /xm surface brightness and that sample dif- 
fer ent shell morpholog ies (disk, shell and torus, 
see lMizuno et al.ll2010l for more details about the 
morphologies). However, we do not consider the 
4 MBs presented in this paper represent to any 
extent a sample from which we could derive con- 
clusions for the whole catalog of 428 MBs detected 
in the Galactic plane. The selected MBs are given 
in Table [T] along with their identification number, 
name, coordinates, size and fluxes. Only MB4121 
is detected in the MIPS 70 /im images and has a 
central source detected in the MIPS 24 /im images. 
Both MB3957 and MB4121 show a central source 



in 2MASS and IRAC images. 

We select the short and long wavelengths high 
resolution {R ^ 600) modules (SH from 9.9 to 
19.6 /j,m, LH from 18.7 to 37.2 /im) in staring 
mode. For each MB, we have an ON position cen- 
tered on the peak of the MIPS 24 ^m extended 
emission and a dedicated OFF position for back- 
ground subtraction that is observed in the same 
sequence as the ON position. For MB4121, we 
have an additional ON position centered on the 
central source detected in the MIPS 24 /im image 
and observed in the same sequence too. For the 
three brightest MB, we integrate for 4 cycles of 
6 seconds in both LH and SH. For the faintest. 
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(a) 



(b) 



(c) 



(d) 

Fig. 2.— MIPS 24 image (2.1'x2.1') with the positions of the IRS sht for (a) MB4001, (b) MB4006, (c) 
MB3957, (d) MB4121. 



Bubble #4001 



ubble #4006 



Bubble #3957 




(d) 

Fig. 3.— MIPS 24 radial profile of (a) MB4001, (b) MB4006, (c) MB3957, (d) MB4121. The black dots 
are all the pixels within 50"of the MB and the red line is the median as a function of the distance. 



MB4001, we integrate for 2 cycles of 30 seconds 
and 4 cycles of 14 seconds in SH and LH respec- 



tively. The position of the SH and LH slits are 
shown on top of MIPS 24 fim images in Figure 
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2(a)|[2(b)| [2(c)] and [2(d)| 

3. Data reduction 

The first step of our data reduction corresponds 
to the co-addition of the different cycles of obser- 
vations (basic cahbrated data or BCD). We use for 
each pixel the median of the stack of BCDs. We 
build a master 2D image for both on target and 
background positions separately. We then sub- 
tract from each master image the corresponding 
master background. This process enables us to 
clean some of the rogue pixels. It also removes 
all the emission features from the Galactic ISM. 
In particular, we detect several weak emission fea- 
tures including II2 lines and PAH features with the 
exact same intensity towards the background and 
the MBs. Consequently, these features which are 
commonly observed in a crowded and dense region 
such as the Galactic plane disappear from the MBs 
spectra after the background subtraction. Other 
gas lines are detected in the master background 
with fluxes at least a 100 times fainter than in the 
master image. We then perform another cleaning 
of the 2D images thanks to the IDL package devel- 
oped by Jim Ingalls (SSC): IRSCLEAInET We use 
SPICEu to extract each background-subtracted 
spectrum. The IRS orders are overlapping each 
other over a few tenth of microns. We use these 
overlaps to manually trim the few first and/or last 
pixels of the order that are usually too noisy, un- 
less both orders appear to be satisfying, in which 
case we average them. 

The slit size of the 2 high resolution modules 
are different (11x22" for LH, 4x11" for SH) but 
both are similar to the characteristic size of our 
targets at 24 /xm. As a consequence, we can as- 
sume our sources to behave neither as point source 
nor as extended source for the spectrum extraction 
purpose in SPICE. Therefore the most adapted 
calibration, which is wavelength dependent, is un- 
certain. We take into account the difference be- 
tween those two extractions when estimating un- 
certainties on gas line fiuxes. There is also a re- 
sulting mismatch (about a factor 2) between the 
SH and LH module spectra, even after correcting 
for the difference in the slit surface. The discrep- 
ancy might be caused by the fact the SH and LH 



slit do not probe exactly the same regions (see 
Fig. [2]). In order to correct these issues, we first 
scale the SH module to the LH module so that the 
dust continuum appears continuous in the overlap- 
ping wavelength range (19.0 to 19.5 /im). We scale 
the entire SH module, i.e. its continuum and the 
gas lines, the same way. We do this only for the 
two MBs that show a significant dust continuum 
in their spectrum. For the two other MBs, that 
exhibit a weak continuum, the match between the 
SH and LH modules might not be correct. While 
ratios between gas lines from the SH and the LH 
module might be inconsistent with models, they 
could be used to constrain the scaling factor be- 
tween the two modules. We discuss this in section 
14.21 We then scale the whole IRS spectrum of 
all four MBs so it matches the background sub- 
tracted MIPS 24 ixm flux as seen through the LH 
slit and as measured in the MIPSGAL images. All 
the scaling factors that we apply are wavelength 
independent and are given in Table [21 The varia- 
tions due to the nodding remain a negligible effect 
(about 1%) except for the outer sheU of MB4121. 
The scaling to MIPS 24 jxm is closer to unity (1.26 
to 1.35) with a point source extraction for the two 
MBs that are the most compact and closer to unity 
(0.87 to 0.91) with an extended source extraction 
for the two MBs that are the most extended. 

4. Results 

The IRS high resolution spectra of the four MBs 
in our sample are shown in figure |4R||4(b)| [8land 
[T3I We classify them in two categories. The first 
corresponds to dust-poor bubbles and is comprised 
of MB4001 and MB4006. The second corresponds 
to dust-rich bubbles and is comprised of MB3957 
and MB4121. The goal of this paper is to identify 
the nature of these objects thanks to the determi- 
nation of some physical characteristics. To do so, 
we use MAPPINGS Ilfl, a code of photoioniza- 
tion, and a simple model for the MBs. 

4.1. The use of MAPPINGS III 

In order to quantify the excitation conditions 
within our MBs, we use the MAPPINGS III code 



^http: / /ssc. spitzer.caltech.edu/dataanalysistools/tools/irsclcan/ 
•^http: / / ssc.spitzer.caltech.edu/dataanalysistools / tools / spice / 



^Modelling And Prediction in Photoloniscd Nebulae and 
Gasdynamical Shocks, developed principally by Ralph 
Sutherland, Mike Dopita and Luc Binette and available 
at |http://www.mso.anu. edu.au/~ralph/map. html| 
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Table 2: Correction factors used to match the SH and LH modules in the case of the 2 MBs that show 
a continuum, and to match the MIPS 24 fim flux for all MBs. The correction factors are defined as: 
Ii, = C'Fi X {Iy{SH)/CF2 + I^{LH)) where 1^ is the final corrected spectrum, and Iy{SH) and Iy{LH) the 
SH and LH spectra extracted in SPICE, either with a point source (ps) or extended source (es) calibration. 
The given uncertainties result from the comparison between nodding positions. 



C Flips) CFijes) CF2{ps) CF-ijes) 



MB4001 


1.35 


± 


0.04 


1.91 


± 


0.05 






MB4006 


1.26 


± 


0.01 


1.77 


± 


0.01 






MB3957 


0.63 


± 


0.02 


0.87 


± 


0.02 


2.06 ±0.04 


1.44 ±0.03 


MB4121 (central source) 


0.65 


± 


0.01 


0.91 


± 


0.01 


2.6 ±0.1 


1.82 ±0.06 


MB4121 (shell) 


0.63 


± 


0.01 


0.87 


± 


0.01 


2.1 ±0.4 


1.5 ±0.3 



of photoionization and model the MBs as iso- 
choric spheres of gas heated by a central blackbody 
source and ionization bounded. We use a template 
provided with the MAPPINGS III code. Details 
about the template are given in Table [31 The MBs 
are modeled as isochoric sphere of gas illuminated 
by a source of 10^ km radius. We use otherwise 
solar abundances for the gas composition, except 
for iron. While iron may not be present in evolved 
low mass stars, it may play an important role in 
the cooling of the more massive ones. For similar 
reason, we also decide to run the MAPPINGS III 
code with and without dust. Therefore we only 
add iron and dust in the model for the two MBs 
that show both iron lines and dust continuum in 
their IRS spectrum (see section 14.3.11 and I4.3.2p . 
The parameters we use for the dust are also given 
in Tabled 

In the following sections, we compare the gas 
line fiuxes predicted by the code to those observed 
in the MBs to constrain, among the many param- 
eters available in the code, the electronic density 
of the bubble and the temperature of the inner 
source. To do so, we run MAPPINGS III with 
density ranging from 100 to 15000 cm^"^ and in- 
ner source temperature from 10** to 3 x IQ^ K. 
The characterisation of the inner source then pro- 
vides constraint on the nature of the MB. We also 
discuss how the presence of dust may affect the 
results. 



4.2. Dust-poor objects: 
MB4006 



MB4001 and 



and 2(b)). Their radial profile at 24 /im exhibit 
only a slight difference at the center: MB4001 is 
flatter and MB4006 has a deeper central hole. For 
this reason they have been classified as a "disk " 

(|2010l l. 



and a "ring" respectively in iMizuno et al. 
Their IR spectrum are also very similar and are 
shown in Figure |4(a)| and [4(b)| They have a very 
weak continuum and many high ionization gas 
lines ([O IV] 25.9 /im, [Ne V] 24.3 and 14.3 ^m, 
[Ne III] 15.5 p.TCi, [Ar V] 13.1 p^va, [S III] 18.7 and 
33.5 /xm). The fiux of these lines are given in Ta- 
ble 2] and [5l We first identify those two MBs as 
planetary nebulae candidates and then constrain 
some of their properties. 

4-. 2.1. Identification as PNe candidates 

These two MBs sha re common properties with 



The first two MBs of our program are morpho- 



logically very similar at 24 /im (see Figure 2(a) 



SMP83 in the LMC (|Bernard-Sala s et al.Jj200J) 
and M l-42 in the Galactic center ( ^ottasch et al.l 
20071 ). two known planetary nebulae (PNe). They 
are also ver y similar to t he ob ject discovered in 
Gepheus by Morris et al.l ( 20061) and suggested as 
a dust-free supernova remnant: detection within 
MIPS 24 fixn bandwidth only (neither in IRAC 
3.6 to 8.0 fim nor in MIPS 70 and 160 ^m), high 
excitation gas lines ([O IV] 25.9 /xm, [Ne V] 24.3 
and 14.3 /im), very little dust continuum. How- 
ever, the gas lines ratios are significantly different. 
In particular, our dust-poor MBs exhibit I([Ne V] 
24.3 fim) > I([Ne III] 15.5 /tm) and I([Ne V] 
14.3 yLtm) > I([Ne III] 15.5 /xm), contrary to the 
Gepheus bubble. 

To accou n t for the lack of dust emission, 
Morris et al. ( 20061 ) suggest that very little in- 
teraction between the shell and the ISM has oc- 
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Fig. 4.— IRS spectra of (a) MB4001 and (b) MB4006 with detected gas lines indicated. 



curred and that the shell has a very high gas-to- 
dust ratio. They rule out AGB, post-AGB and 
PN shell based on the following arguments: (1) 
the electronic density they derive from [S III] 
33.5 /im to [S III] 18.7 fim does not depend on 
the radius while multilayered shell would create 
large variations in N^, (2) typical IR spectra of 
PNe show PAHs features and higher ionized gas 
species, (3) the nebular gas spectrum should con- 
tain hydrogen lines tracing the hydrogen lifted 
away from the surface of the star in its AGB 
or post-AGB phases prior to the star emerging 
the H-deficient white dwar f cooling track. How- 
ever, iFesen fc Mihsavlievid (|2010l ) obtained nar- 
row and broad passband optical images as well as 
low dispersion optical spectra of the Cepheus bub- 
ble. They reject the dust-free SNR interpretation 
based on the absence of high- velocity features, the 
presence of relatively strong [N II] emission and 
the lack of detected [S II] emission, which would 
indicate the presence of shock-heated gas. Addi- 
tionally, no coincident X-ray or nonthermal radio 
emission has been detected. They suggest it is 
more likely a faint high-excitation PN, lying at a 
distance ^ 2.5 kpc with strong [O IV] emission 
rather than dust continuum emission dominating 
in the MIPS 24 /im images. They have identified 
a possible central star candidate with m^' ~ 22.4. 

Furthermore, while the presence of PAHs 
features in PNe has be en reported by man y 
authors: in NGC7027 dBeintema et al.l Il996l) . 



Magellanic Clouds PNe (several exampl es in 



Stanghellini et al.ll2007l:lBernard-Salas et al.ll2008l 



20091 ) and Galact ic Bulge PNe (one ex ample. 



PNG351.2-h05.2 in iGutenkunst et al.ll2008l) . there 
are also examples of PAH free PNe, sometimes 
in the same references: S MP83 in the Large 
Mage llanic Cloud (LMC) ()Bernard-Salas et al. 
20041 ) ■ a rather faint PN in the direction of the 



20071) ■ 



Galactic center (Ml-42 IPotta sch et al 
Magellanic Clouds PNe (several exampl es m 



Stanghellini et al.ll2007l:lBernard-Salas et al.ll2008l 



200l and Galactic B ulge PNe (ah but one in 
iGutenkunst et al.ll2008i ). Among those PNe that 
show PAH-free IR spectrum, some present a signif- 
icant continuum at wavelengths longer than 15 jum 
or crystalline silicate feature s ( Stanghellini et all 
20071 : IGutenkunst et al.ll2008t ) while others hardly 



show any continuum or broad features from 10 to 
35 M m. MB4001 and MB 4006 as well as the bub- 
ble of iMorris et al. (2006) belong to this last class 
of object, which may be called dust-poor PNe. 
Other examples are S MP13, SMP35, SMP40 and 



SMP83 in the LMC ( Bernard-Salas et al.l 12004 , 



2008") as well as Ml-42 toward the Galactic center 



(iPottasch et al.ll2007l) 



Finally, iGroves et all (|2008l ) studied the neon 
and sulfur emission lines for a wide variety of 
objects (e.g. galaxies, PNe, HII regions) from 
the Spitzer and ISO archives. Their diagram of 
[Ne III] 15.6 fim /[Ne II] 12.8 fim vs [S IV] 10.5 
/[S III] 18.7 /im shows a tight correlation. The two 
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dust-poor MBs would fall on the high end of this 
correlation, where the PNe, which represent the 
most isolated population in their sample, are the 
most numerous. 

We conclude that MB4001 and MB4006 are un- 
doubtedly PNe candidates. 

4-. 2. 2. Characterization 

Even though the examples of dust-poor PNe 
in the literature are not numerous, MB4001 
and M B4006 both show strong s imilarities with 
SMP83 (|Berna rd-Salas et al."2004') and, to a lesser 
extend, Ml-42 (P ottasch ct al. 2007) and Sh 2-188 
()Chu et al.ll20C)i . We compare gas lines fluxes of 
the two MBs, normalized to [S IV] 10.5 /xm in the 
SH module and to [S III]33.5 /im in LH module, to 
those of SMP83 in Figure[5l We normalize and dis- 
cuss both modules separately, as some mismatch 
between the modules might remain (see section 
[3]). Within the SH module, we observe that the 
[S IV] 10.5 /im to [S III] 18.7 /im line ratio appears 
to be very similar within the three objects. The 
same can be told about the [Ar V] 13.1 /xm line 
relative to the [S IV] 10.5 /xm line. The other lines 
within the SH modules are three Neon lines at 
three different ionization levels ([Ne II] 12.8 /im, 
[Ne III] 15.5 /xm and [Ne V] 14.3 ^.m). Based on 
their relative fluxes MB4001 and MB4006 appear 
to be more highly ionized than SMP83. Within 
the LH module, all the gas lines arise from differ- 
ent elements. On the one hand, the ratio between 
[Ne V] 24.3 /im and [S III]33.5 /xm is consistent 
within the three objects. On the other hand, 
the [Si II]34.8 /xm line is much weaker, relative 
to [Ne V] 24.3 /xm, in MB40Q6 than in SMP83. 
It is not detected in MB4Q01. Additionally, the 
[O IV] 25.9 /xm appears to be stronger in the 
MBs than in SMP83. We note that it's again 
within MB4Q01 that excitation, as traced by the 
[O IV] 25.9 /im line, seems to be the strongest. 

In order to quantify the excitation conditions 
within the two dust-poor MBs, we use MAP- 
PINGS HI with the parameters presented in Ta- 
ble |3] and without dust. We use the results 
from these models to compare the IR gas line 
fluxes to those observed in the dust-poor MBs 
and then constrain the temperature of the inner 
source of radiation and the gas density of the 
MB. In order to avoid the effect of any remain- 
ing LH/SH mismatch and any dependence on the 
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Fig. 5. — Comparison of gas line fluxes between 
MB4001 (red), MB4 06 (g reen) and SMP83 from 
Bernard-Salas et al.l (2004) in diamonds. Fluxes 
are normalized to [S IV] 10.5 /xm = 1 in SH module 
(left of the dashed line) and to [S III] 33.5 /xm =0.1 
in LH module (right of the dashed line). 



chemical composition of the gas, we use as diag- 
nostics the [S IV] 10.5 /xm to [S III] 18.7 /xm and 
[Ne V] 14.3 /im to [Ne III] 15.5 /xm hue ratios. Fig- 
ure [6] shows the evolution of both these ratios as 
predicted by our models for various temperature 
and density. Measurements for MB4001, MB4006 
and the planetary nebula SMP83 are also shown. 
The line ratios of SMP83 are well reproduced by a 
density of 5000 ± 2500 cm^"^ and a source temper- 
ature of about 170000 ± 10000 K. The density we 
estimat e is in agreement with pr evious measure- 



ments ( Bernard-Salas et al.l 12004 and references 



therein). The temperature of the c entral source 
i s the same as that estimated by iDopita et al 



(|1993[ ) with MAPPINGS 2 and a multizone model 
combining an optically thick ring with optically 
thin extensions. They found densities of 1000 and 
3600 cm~^ in the thin and thick components re- 
spectively. The required densities are 10000 ± 
5000 cm-3 and 5000 ± 1500 cm'^ in MB4001 and 
MB4006 respectively, which are similar to that of 
SMP83. The inner source temperature is signifi- 
cantly hotter though. The IR gas lines of MB4001 
are in agreement with T = 205000 ± 5000 K 
while those of MB4006 are in agreement with 
T = 225000 ± 10000 K. Those temperature are 
on the high-end of what is found in published ex- 
amples of PNe. 
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Using these parameters, we are able to predict 
ratios between SH and LH lines (e.g. [S III] 33.5 /xm 
to [S III] 18.7 and [Ne V] 14.3 /zm to 

[Ne V] 24.3 fim) and thus, by comparing them 
to the observed fluxes, estimate the scaling factor 
that should be applied to match the SH module 
with the LH module. For MB4001, the low end 
of the density range (~ 5000 cm~'^) leads to an 
SH scaling factor of about 4.5 and 7.0 using the 
[S III] and [Ne V] line ratios respectively, while 
the high end of the density range (~ 15000 cm~'^) 
leads to a scaling factor of about 8.4 and 17.0 us- 
ing the [S III] and [Ne V] line ratios respectively. 
For MB4006, those values are 3.7 and 5.5 at low 
densities 3500 cm~'^) and 6.3 and 15.0 at high 
densities 6500 cm"'^). We decide not to apply 
any scaling factor to the SH module given such a 
large uncertainty. However, we draw the reader's 
attention to the underestimation, by a factor of a 
few, of the SH lines fluxes given in Table U and [S] 

From the comparison of the gas line fluxes with 
those of SMP83, we also give a rough estimate 
of the distances to MB4001 and MB4006. We 
use only LH gas lines to perform this comparison. 
The [Ne V] 24.3 and [O IV] 25.89 fj.m are respec- 
tively 4.3 and 14 times brighter in MB4001 than in 
SMP83. Those values are 23 and 28 for MB4006. 
The distance to SMP83 is 50 kpc. To scale the 
emission lines of MB4001 and MB4006 adequately 
requires them to be between 18 and 200 times 
closer (0.25 and 2.8 kpc) and between 530 and 
780 times closer (0.06 and 0.1 kpc) respectively, 
depending on the gas line used for the scaling. 
These estimations assume the three PNe have sim- 
ilar physical properties. However, since we know 
the two MBs are significantly hotter than SMP83, 
those distances might be lower-limits. Another es- 
timate can be provided by the apparent sizes of the 
MBs. At a distance of 50 kpc, SMP83 is smaller 
than 1", at least 15 times smaller than MB4001 
and 21 times smaller than MB4006. Therefore, 
the two MBs would be at a distance of at most 
3.3 and 2.4 kpc respectively. At those distances, 
the MBs would have a physical size of about 0.1 to 
0.2 pc which is satisfying. Moreover, at an arbi- 
trary distance of 2 kpc, a white dwarf of 200000 K 
surface temperature and 5000 km in radius has an 
apparent magnitude of at least 20 from V band to 
IRAC and MIPS channels. Therefore such tem- 
peratures and distances are consistent with the 



non-detection of the inner source of the dust-poor 
MBs, especially in the Galactic plane where the 
extinction is high. 

Finally, we discuss the very limited detection 
of dust emission in the MBs spectra. No PAH 
features are detected at short wavelengths in the 
IRS spectrum and the IRAC images. The MIPS 
70 /im broadband images (not shown here) do not 
present any trace of colder dust. The MIPS 24 /im 
is mainly accounted for by the [Ne V] 24.3 /xm and 
[O IV] 25.9 ^m hues. In the case of MB4001, a 
weak continuum (< 0.1 Jy) is detected in the LH 
range of the spectrum but its contribution to the 
total MIPS 24 fim flux is only about 30%. In the 
case of MB4006, the continuum is hardly detected 
and its contribution to the MIPS 24 /im emission is 
less than 15%. The detection of a weak [Si II] line 
in MB4006 indicates that a faint PDR might still 
be present in this object. As we already mentioned 
in section 14.2.11 examples of such dust-poor PNe 
exist in the literarure. The most likely explanation 
is that the dust and the PDR have been destroyed 
by the hard radiation field of the inner source. It 
has been found that larger, more evolved, PNe 
usually exhibit less dust continuum and more 
[O IV] gas line ernission in their MIR sp ectrum 

' 20091) as the 



( Stanghellini et al.ll2007t iChu et al 



amount of ionized gas increases at the expense of 
the molecular and atomic material ("Hug gins et al.l 
1996; Bcrnard-Salas & Ticlens 2005). In this 
global picture, MB4001 and MB4006 would be 
template of old PNe. 

4.3. Dust-rich objects 

The other two objects in our sample share at 
least two properties: (1) a IR central source is 
detected in several Spitzer channels and (2) their 
mid-IR spectra is dominated by a continuum. 
However, these two objects are as much different 
from the first two as they are from each other, 
both in terms of IR morphology and spectroscopy. 
We thus present them separately. 

4.3.1. MB3957 

The MIPS 24 /xm morphology of MB3957 is 
that of a torus. A lack of MIR emission at the 
center is obvious from the MIPS 24 ^m image 
and radial profile (see Fig. 2(c) and 3(c) I. Fig- 
ure [7] shows the IRAC observations of MB3957, 
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(Ne v)14.3 / (Ne lll]15.5 



Fig. 6.— [S IV] 10.5 fiui to [S III] 18.7 fim ratio as 
a function of [Ne V] 14.3 /xm to [Ne III] 15.5 
ratio from the MAPPINGSIII model without dust 
for different set of temperature of the inner source 
and density of the isochoric sphere of gas. From 
bottom to top, the curves are representing den- 
sities of 1000, 3000, 5000, 7000, 9000 and 10000 
cm~'^. From left to right, diamonds are for tem- 
peratures of 150000, 200000, 250000 and 300000 
K. MB4001 and MB4006 ( this paper) and SMP83 



( Bernard-Salas et al.ll2004D are indicated with the 



red, green and blue crosses respectively. 



from the GLIMPSE survey (|Beniamin et al.ll2003[ ) 
with contours from the MIPS 24 /xm observation. 
Extended emission associated with MB3957 is vis- 
ible at every IRAC wavelengths, though the pres- 
ence of nearby bright diffuse emission prevents us 
from an unequivocal conclusion. A point source 
is detected at the center of MB3957 in all IRAC 
bands as well as 2MASS J, H and K bands (not 
shown here) but not at 24 /im. It is identi- 
fied as "SSTGLMC G305.6517+00.3493" in the 
GLIMPSE I Spring '07 Catalog and as "13145704- 
6223533" in the 2MASS All-Sky Point Source Cat- 
alog. Based on the 2MASS and IRAC colors of 



the central source. iMarston et a l. (2010) have sug- 
gested MB3957 as a Wolf-Rayet candidate. The 
ON position of the IRS slit only partially covers 
the central source in both the SH and LH modules 



(see Fig. [2(c) I. 

The Spitzcr/IRS spectrum of MB3957 is pre- 
sented in Figure [8l It is dominated by a contin- 
uum and several low excitation gas lines whose 
fluxes are given in Table |6l At short wavelengths, 






Fig. 7.— IRAC 3.6 (top-left), 4.5 (top-right), 5.8 
(bottom-left) and 8.0 /im (bottom-right) images 
(2'x2') with contours of MIPS 24 ^m of MB3957. 

no PAH features are visible. As shown on the 



MIPS 24 ^m and IRAC observations (see Fig. 2(c) 
and [7]) , the background emission is varying sig- 
nificantly around MB3957 and the IRS OFF po- 
sition is in the lowest surface brightness region. 
We check that the background surface brightness 
is not varying enough in the MIPS 24 and IRAC 
8 /im images to account for the residual contin- 
uum in the final IRS spectrum of MB3957 and 
that the MIR continuum is actually arising from 
the MB. In the following, we combine the infor- 
mations at our disposal (central source photome- 
try and MIR spectroscopy) to constrain the true 
nature of MB3957. 

4.3.1.1. The MIR gas lines The MIR spec- 
trum of MB3957 shows several gas lines listed in 
Table [6] These lines are not resolved with IRS 
(i? ~ 600). We correct the line fluxes for extinc- 
tion. The details of the correction are given in 
section liXOl The detection of [Fe III], [Ne II], 
[Ne III], [S III] and [S IV] hues, along with the non- 
detection of, for instance, [Fe II] , [O IV] or [Ne V] 
lines, sets some constraints on the energy of inci- 
dent photons. Table [7] gives the ionization poten- 
tials of Iron, Neon, Sulfur and Oxygen for their 
ionization states that emit in the MIR. There is a 
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Fig. 8.— IRS spectrum of MB3957. 

range of energy where these elements with these 
observed ionization states can coexist. Photons 
with energy ranging from 16 to 41 eV have to be 
significantly more numerous than those with en- 
ergy below 16 eV or above 55 eV. Simply based 
on these ionization potentials, we expect the shell 
to emit [O III] 51.8 and 88.4 fim lines. 

We first discuss the detection of the [Fe III] 22.92 fim 
line. We start by verifying that it is not the [Fe II] 
22.90 fiia line. The previous discussion on the ion- 
ization potentials already suggests that the pres- 
ence of Fe III is more likely than that of Fe II. 
Then, if it was the [Fe II] 22.90 /im line, our mod- 
els predict we should also see the [Fe II] 17.9 and 
26.0 iim lines with fluxes at least an order of mag- 
nitude higher. Since we do not detect those lines, 
we conclude that the [Fe II] 22.90 /im line is at 
least a factor 10 fainter than the sensitivity in the 
LH module which we estimate to be no better than 
4 X 10~^^ W.crn"^. As a consequence, the [Fe II] 
22.90 fim line contributes at most 1% to the fea- 
ture at 22.9 nm. Therefore, it is more likely that 
this line is the [Fe III]] 22.92 ^m hue. This line 
has rarely been observed: in the supernova rem- 
nant RCW103, prototype of a SNR shock hea vily 
interacting with dense ISM ( Oliva et"al] 19981). in 



the PNe NGC7009 and NGC7027 (Rubin et al 



1997), and toward the Galactic center (jContini 
20091 ). [Fe III] lines are more commonly observed 



the [Fe II] lines at 17.9, 24.5, and 26.0 /xm are usu- 
ally significantly brighter than the [Fe III] 22.9 ^m 
line. In MB3957, the [Fe III] 22.92 /xm line is at 
least a factor 10 brighter than the [Fe II] 24.5 and 
26.0 /im lines. Indeed, we estimate the sensitivity 
to be no better than 4 x 10~^^ W.cm^^ within the 
LH module, where the [Fe II] 24.5/26 fim, [O IV] 
25.9 /im and [Ne V] 24.3 //m would be if they were 
detected. 

Then, the detection of an iron line might trace 
the presence of a shock that processes the dust 
grains and sends iron atoms back to the gas. In 
the case of a shock, the sulfur and neon gas lines 
are in agreem ent with an electronic density of a 
few 10^ cm~'^ ( Alexander et al.| 1999 ) and a veloc- 
ity of 200 - 300 km.s"^ (|Hewitt et al.ll2009f) . The 
MAPPINGS I II Library o f Fast Radiative Shock 



in the optical part o f the spectrum, even though 
Keenan et all (|l992[ ) showed the [Fe III] 22.92 /im 
was the strongest line of the Fe 3d^ configuration 
in a 10^ K and 100 to 1000 cm'^ gas. Moreover, 



Models from lAllen et all ([2008!) also predicts 
that, within a density of 100 cm~'^, the observed 
[Ne III] 15.5/[Ne II] 12.8, [S III] 18.7/[S III]33.5 
and [S IV] 10.5/[S III]33.5 are in agreement with a 
shock of 200 — 250 km.s^^, along with its precur- 
sor. However, one would expect to see other MIR 
iron lines (e.g. [Fe II] 17.9, 24.5 and 26.0 /im). In- 
deed, the fast radiative shock models also predict 
that the [Fe II] lines are at least one order of mag- 
nitude brighter than [Fe III] 22.92 /im. Therefore 
we suggest that photoionization rather than shock 
is the interpretation for the neon, sulfur and iron 
lines. 

We then use the grid of models that we com- 
puted with MAPPINGS III. Figure [9] shows two 
diagrams of the [S III] 18.7/[S III]33.5 ratio as a 
function of the [Ne II] 12.8/ [Ne III] 15.5 ratio for a 
given inner source temperature or a given gas den- 
sity. The line ratios observed in MB3957 match a 
central source temperature of 60,000 to 70,000 K, 
with a gas density of about 250 to 3000 cm~'^. 
These results are not significantly altered by the 
addition of dust in the model. While significantly 
cooler than the central source of the two dust-poor 
MBs, the source at the center of MB3957 is sig- 
nificantly hotter than a main sequence star and 
would match a white dwarf or a Wolf-Rayet star. 
MAPPINGS HI does not include the [Fe III] lines 
fluxes, but it gives the relative abundances and 
total column densities of the iron first ionization 
levels. In order to have more Fe III cations than 
Fe II and Fe IV cations, the temperature of the 
central source has to be between 30,000 and 80,000 
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Fig. 9.— Diagrams of the [S III] 18.7/ [S III]33.5 
ratio as a function of the [Ne II] 12.8/[Ne III] 15.5 
from our MAPPINGS III modehng. (a) shows 
curves for a given gas density from 250 to 15000 
cm~'^. (b) shows curves for a given inner source 
temperature, from 50000 to 150000 K. The black 
cross shows the observed ratios for MB3957. 



K, which is in agreement with the neon and sulfur 
lines diagnostics while the gas density is not well 
constrained. The maximum Fe III-to-Fe II abun- 
dance and column density ratios are both about 
3.4 for iTBB,ngas) = (60000 if, 3000 cm^^) and 
(50000 if, 500 cm^'^) respectively, a factor of at 
least a few smaller than the observed flux ratio be- 
tween the [Fe III] 22.9 ^m line and the [Fe II] 24.5 
and 26.0 ^m lines. To explain the remaining dis- 
crepancy, we could invoke the difference between 
the oscillator strength of the [Fe III] 22.9 fim line 



and that of the [Fe II] 24.5 and 26.0 /zm lines. Un- 
fortunately, these values are not available due to 
the complexity of the calculations for the Fe III 
cation. Difference of extinction between the iron 
cations emitters, due to layering of photo-ionized 
matter in the shell could be another interpretation 
of the remaining discrepancy since Fe II cations 
are expected to be farther away from the central 
source and deeper in the shell structure than the 
Fe III cations. However, to attenuate the [Fe II] 
lines by a factor 3 larger than the [Fe III] lines re- 
quires an extinction differential of about 1.1 mag- 
nitude at wavelengths of about 20 /xm, which cor- 
responds to about 40 magnitudes of visual ex- 
tinction. This seems unrealistic. The presence 
of the [Fe III] 22.9 iim hne without the mid-IR 
[Fe II] lines in the spectrum of MB3957's shell 
remains complicated to interpret due to the dif- 
ficulty to model with precision the Fe III cation 
transitions. However, we suggest the whole set 
of observed mid-IR gas lines in MB3957 indicate 
that the shell is produced by the photo-ionization 
of a few 10^ cm~^ gas by an inner source of about 
6 X 10'' K. 

4.3.1.2. The central source The 2MASS 
and IRAC fluxes of the central source in MB3957 
are given in Table [8] and are shown in Figure 1101 
The modehng of the MIR gas hues with MAP- 
PINGS HI suggests that the central source has 
a temperature of about 60000 K, which could fit 
either a white dwarf or a very hot Wolf-Rayet 
star, as suggested by iMarston et al.l (|2010i) . The 
SED peaks in the near-IR which clearly indicates 
a reddened source. Indeed, the best fit of the en- 
tire SED with a blackbody leads to a temperature 
of 1600K (see Fig. 10(a)). Reddening might be 



due to extinction along the line of sight, highly 
plausible through the Galactic plane. However, 
with an intrinsic temperature of 60000 K for the 
source, it's impossible to find a decent fit of the 
SED, whatever the distance and the extinction. 
We thus suggest the IR excess in the SED orig- 
inates in a circumstellar enveloppe of dust or in 
the presence of a substellar companion around the 
inner source. The presence of a substellar com- 
panion in a binary system t hat comprises a white 
dwarf has been reported bv lZuckerman fc Becklin 
(|l992t ). GD165 and GD1400 are two exam- 
ples of such systems (jKirkpatrick et al. I1999I: 
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Farihi fc Christopheij 120041) . The presence of a 
disk of du st around a wh ite dwarf has been ob- 
served by ISu et all (|2007l ) in the MIPS 24 
image of the Helix Nebula. Wolf-Rayet stars have 
long been known to be surrounded by dust shells. 
We use the SED of the central source in MB3957 
to determine the validity of either interpretation. 

We model the central hot source and the cooler 
source (surrounding disk of dust or brown dwarf 
companion) with blackbody functions. The best 
fit depends on the initial guess and usually leads to 



a one component fit (see Fig. 10(a) I. The reduced 



for this fit is about 9. In order to "force" the 
use of two components, we fix the dilution factor 

to non-zero 



of the central source 



values and run several fits. The other parameters 
(star temperature, dust temperature, dust scal- 
ing factor, extinction scaling factor) are set free 
with some boundaries: 50000 K < T* < 70000 K, 
10 K < Tdust < 10000 K and 1 mag < Ay < 
100 mag. We use the average NIR ex tinction 
values found by llndebetouw et al. (l2005l) in the 
Galactic plane. We explore a very wide range of 
values for (i?*/d*)^. Small values of ('^ 10) 
are found for very small values of (i?*/fi*)^, which 
corresponds to the fit with only one component. 
However, we find a better fit {x^ ^ 2) for a very 
narrow range of (i?*/d*)^ (2 — 3 x 10^^^). The cool 
component temperature and the extinction along 
the line of sight both reach well determined values 
in that range: Tdust 900 K, Ay 20 mag. The 
best fit is shown in Figure [10(b)] The dilution fac- 
tor for the best fit (i?,/d* = (2.5 ± 0.5) x 10"") 
allows us to constrain the nature of the central 
source. If the central source is a white dwarf 
(i?* ~ O.Oli?©) it is at a distance of 4.5 ± 0.5 pc 
which is highly unlikely. Indeed, the source would 
be one of the closest star to Earth. Moreover, 
the best fit extinction thus corresponds to a col- 
umn density of 3.7 x 10 ^^ cm~^, using Ay/Nn 



0.53 x 10-21 cm-2 from lSavage fc Mathij (Il979l) . 
Therefore, the average density along the line of 
sight would be about 3 x 10"^ cm"'^, several order of 
magnitude larger than that of the diffuse medium. 
If the central source is a WR star (i?* ~ lOi?©) 
it is at a distance of 4.5 ± 0.5 kpc which is much 
more likely as it corresponds to Galactic distances. 
Moreover, the average density would be about 
3 cm""^ which is very similar to that of the Galac- 
tic diffuse ISM. We thus conclude that the central 



source of MB3957 is a WR star and that the red- 
dening of its SED is due to the presence of an in- 
ner circumstellar enveloppe of hot dust. Our con- 
clusion therefo r e confi rms the interpretation from 
Marston et al. ( 2010t ). We characterize the dust 



component in the following section. 

4.3.1.3. Full SED fitting The complete SED 
of MB3957, combining the 2MASS-hlRAC SED 
and the IRS spectrum, is shown in Figure [TT] with 
an intensity scale that allows a better view on the 
continuum emission. The 2MASS and IRAC ob- 
servations are well interpreted by the two compo- 
nent fitting presented in the previous section. In 
the IRS wavelength range, the emission is domi- 
nated by a continuum that shows silicate absorp- 
tion features at 10 and 20 /zm and that dominates 
the emission in the MIPS 24 fj,m. Therefore, we 
believe the entire SED of MB3957 can be inter- 
preted with three emission components: the cen- 
tral source, an inner disk of hot dust and an outer 
shell of warm dust. 

In order to confirm this interpretation, we 
model the entire IR SED of the shell and the 
central source with three blackbodies attenuated 
by extinction along the line of sight. We com- 
bine the so-called "ISM " and "Galactic c enter" 
extinction curves from Chiar fc Tielens (|200d) 
into several extinction curves, defined by the rel- 
ative contribution of the "ISM" extinction curve 
Uism)- We set free seven parameters: the tem- 
peratures {T^,,Thot,Tyjarm) and dilution factors 
{W,^,Whot,Wtnarm) of the blackbodics, and the 
normalization of the extinction curve at 9.7 /im 
(rg.y). We put limits on the three temperatures: 
50000 to 70000 K for the central source, 500 to 
2000 K for the hot dust and 10 to 200 K for the 
warm dust. We then run a grid of fits with the pa- 
rameters W* and fisM fixed at pre-defined values: 
fiSM varies between and 1 and varies from 
IQ-^^ to 10-2°, following the results of previous 
section. The lower values of are obtained for 
fiSM ~ and for W^, ^ 0.9 x lO^^i^ which is only 
a factor of a few smaller than the value found for 
the fit of the 2MASS and IRAC SED only. As 
a consequence, the conclusion on the nature of 
the central source remains the same: MB3957 is a 
WR. The distance to MB3957 is 7.5 x (R^/IORq) 
kpc. We do not show here the contour plot for the 
other parameters but we discuss their values and 
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Fig. 10. — One and two component fitting of MB3957 central source SED. The black crosses are the fluxes 
from the 2MASS and GLIMPSE catalogs, (a) Best fit with a single component (red diamonds), (b) Best- fit 
of the central source SED (red diamonds show the cold component, blue diamonds show the hot central 
source and black diamonds show the total) with two components. 



variations. For the best fit, the depth at 9.7 ^m 
and the temperatures of the hot and warm com- 
ponent are 2.0, 1750 K and 167 K respectively. 
The warm component temperature does not vary 
significantly within the grid we explore, from 165 
to 185 K, as it is well constrained by the IRS spec- 
trum. The hot component temperature increases 
from 500 to 2000 K as VF* decreases but is no- 
ticeably above 1250 K for W^* < 2 x 10"^^ Those 
values are in agree ment with previous observa tions 



of WR stars (e.g. Ivan der Hucht et al.l 119961 ). At 



those temperatures, silicates would begin to be 
detected in emission rather than in absorption ex- 
cept in the case of carbon rich environment, from 
which silicates would be absent. The depth T9.7 
decreases from 4 to 1.5 as decreases but is 
noticeably below 2.5 for < 2 x 10~^^. Since 
fiSM ^ for the best fit, we use the "Galac- 



tic cen ter" extinction curve from lChiar fc Tielens 
(|2006l) to convert Tg.7 = 2.0 to Ay = 13.7 mag, 
which is in good agreement with the extinction 
obtained in the previous section. The tempera- 
ture of the central source spans the entire range of 
allowed values, varies quickly from 50000 to 70000 
K as decreases from 10"^^ to 3 x 10~^^, and 
for the best fit = 51000 K. All parameters are 
significantly more dependent on than fisM- 

The best fit is shown in Figure [TT] along with 
the dereddened spectrum of MB3957. The warm 



component is by far the main contributor of the 
emission within the IRS wavelength range while 
the NIR data points can almost only be accounted 
for by the central source and the hot component. 
This explains why the results from the fit of the 
IRAC and 2MASS data points only are not signifi- 
cantly different. The shape of the continuum in the 
IRS range is almost perfectly reproduced by the fit 
even though one may argue the silicate feature at 
20 /im seems deeper in the observations. How- 
ever, taking into account the noise in the IRS/SH 
data, and the possible variations of the silicate fea- 
tures in the e xtinction curve not acco unted for by 
the curves of Chiar fc Tielend ()2006h . we believe 
the three blackbody components model is a sim- 
ple but efficient way to interpret the continuum 
emission of MB3957. The dereddened curve, also 
shown in Figure [TTl presents the intrinsic emission 
spectrum of MB3957 without correcting for the di- 
lution factor and scaled up by a factor 10 for more 
clarity. The discrepancy at the very short end of 
the IRS wavelength range might also be explained 
by a lower S/N ratio combined with a higher ex- 
tinction correction factor. 

4.3.1.4. Dust mass loss Here we finally ver- 
ify that the distance to MB3957, assumed to be 
7.5 kpc, its MIPS 24 /im brightness and the fit 
of its whole SED are consistent with the WR in- 
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Fig. 11.— Three blackbody component fitting of the whole SED of MB3957.MB3957's whole SED (black 
line and diamonds), best fit (red lines and diamonds), dereddened SED multiplied by a factor 10 (blue lines 
and crosses) and small dust grains model multiplied by a factor 10 (solid green line). 



terpretation. We model the warm component, 
which accounts for the emission of dust in the 
MIPS 24 /im sheU of MB3957, with a more re- 
alisti c emissivity functi o n. We use the dust model 



from ICompiegne et al.l (|2011l ). For the incident 
radiation field, we use a 50000 K blackbody at a 
distance of 12" to the peak emission of the shell, as 
given by the MIPS 24 /im image and which corre- 
sponds to 0.4 pc at a distance of 7.5 kpc from the 
Sun. We plot the emission spectrum of the only 
small amorphous carbon (SamC) grains in Figure 
[TTl scaled to match the IRS spectrum. The re- 
quired dust column density is 2.7 x 10~^ g.crn"^ 
if we use SamC only. For large amorphous carbon 
(LamC) or amorphous silicate (aSil) grains, the re- 
quired dust column density is'^1.3xl0~^ g.cm~^ 



but the spectral shape of the IRS spectrum is not 
reproduced as well as with SamC. This might be 
an indication that large dust grains are ground 
into smaller one, though a more thorough model- 
ing of the dust size distribution is required. In- 
tegrated over the MB3957, the total mass of the 
shell is 10-3 X [0.01/ Xd) x (d/7.5kpc)2 Mq where 
Xd is the dust-to-gas mass ratio. Therefore, the 
MIPS 24 fim shell of MB3957 would have been 
formed int = 2125xd*(kpc)/(uu,md/500 km.s"^)) 
years where Vwind is the average stellar wind speed 
that produced the shell of the WR star. We use a 
reference velocity of 500 km.s"^ since the mid-IR 
lines are not resolved with the resolution of IRS 
{R ~ 600). The formation timescale is about 2 or- 
ders of magnitude shorter than the mean lifetime 
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of a WR (5 X 10^ yr, iMaeder fc Mevned Il987l) . 



which ahows for muUiple mass foss events. The 
MIPS 24 shell thus traces a very specific time 
in the WR's life. We then infer the average mass 
loss rate of MB3957 to be 4.8 x IQ-^ x (0.01/xd) x 
(d/7.5kpc) X (wu)md/500 km.s^^) Mq/jt, a factor 
of a few lower than the 10~^ to 10~^ M0.yr~^ that 
WR stars have long been known to eject through 
powerful winds (e.g., WiUis 19911). However, the 
shell mass, timescale and wind velocity we obtain 
here are esti mates within a facto r of a few at least. 
For instance, ISahai et al. I (|2nnn[ ) have shown that 
the dust-to-gas mass ratio in the planetary nebula 
with WR type central star He 2-113 is only 0.5%. 
Typical values of 0.25% to 1% are found in the 
circumstellar envelopes of AGB stars, which leads 
to a factor 4 uncertainty. Moreover, we assume 
the hot and warm component are blackbodies to 
estimate the distance to MB3957. If we model 
both with a i/^B^, function instead, the best fit 
(not shown here) is found for slightly different pa- 
rameters. The best fit dilution factor of the star, 
for instance, corresponds to 11 x 10~^^. It still 
implies that MB3957 is a WR star and places it 
at a distance of 2.1 x (R^/IORq) kpc, more than 
a factor 3 closer than the 7.5 kpc determined from 
the three blackbodies fitting. 

Combining all the observational constraints on 
MB3957, from near- to mid-IR broadband pho- 
tometry and mid-IR spectroscopy, we conclude 
that the most likely interpretation for this shell 
is that it is a WR surrounded by two circumstel- 
lar envelopes. The inner one, visible as an excess 
in the IRAC bands, is not resolved and may be 
related to a recent mass loss. The larger one, vis- 
ible in the MIPS 24 fim image, is related to an 
older mass loss. The modeling of the outer shell 
torus-like morphology which has been related to 
the presence of a compa nion in binary systems 
(|Morrialll98it lLiviolll993h . is beyond the scope of 
this paper and would benefit from better angular 
resolution. 

4.3.2. Dust-rich object MB4121 

The morphology of MB4121 is that of a shell 
with a central source in the MIPS 24 /j,m image 
(see Fig. 2(d)). At shorter wavelengths (IRAC 
and 2MASS observations) , only the central source 
is visible (see Fig. [T2|). What looks like diffuse 
emission on the IRAC 5.8 /im is related to instru- 



mental artifact due to the bright central source. 
A nearby IR dark cloud (IRDC) is visible in the 
IRAC and MIPS 24 ^m bands. The central source 
is identified as "SSTGLMC G337.5543-I-00.2198" 
in the GLIMPSE I Spring '07 Catalog and as 
"16364278-4656207" in the 2MASS All-Sky Point 
Source Catalog. The IR magnitude s and fluxes for 
the st ar are reported in Table El IWachter et al 



(l2010[ ) lave recently identified the central source 
of MB4121 as a Be/B[e]/LBV thanks to the strik- 
ing similarity of its near-IR spectrum with that of 
G024.73+00.69. 

Contrary to the previous three objects pre- 
sented in this paper, we observe MB4121 toward 
two positions. The first position is on the central 
source of the shell and the second is on the north- 



ern rim of the shell (see Fig. 2(d) I. The two back- 



ground subtracted IRS spectra are presented in 
Figure [131 A long wavelength (A > 20/im) contin- 
uum is detected toward both positions. Therefore 
we attribute this emission component to the outer 
shell. The spectrum toward the central source 
shows a richer spectrum with a short wavelength 
(A < 20/im) continuum, absorption features of sil- 
icates and CO2 ,and several emission gas lines, 
most of them from iron (see Fig. [13(b)]). The mea- 
sured fiuxes of the main gas lines labeled in the 
IRS spectrum are given in Table [TOl and im 

We first correct for extinction the IRS spec- 
tra using the whole SED of the central source in 
MB4121. We then use the MIR gas lines as a diag- 
nostic to constrain the temperature of the central 
source and the density in the object. We finally 
discuss the outer shell continuum emission. 

4.3.2.1. Extinction correction The spec- 
trum toward the central source shows absorption 
features from silicates at 9 and 18 /im and CO2 at 



15 /im (see figure 13(a) ). These features may arise 



from the IRDC visible in figures 1(d) and [121 at a 



position close to that of MB4121. The outer shell 
of MB4121 may also contribute to the extinction 
along the line of sight toward the central source. 
We correct the emission spectrum of MB4121 to 
take into account the extinction. 

First, we extract the profile of the CO2 bending 
mode at 15 /im foll o wing the method presented in 
Pontoppidan et al. I (I2OO8I) . The continuum is fit- 
ted by a third order polynomial function in the 
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Fig. 13. — (a) IRS spectrum toward the central source of MB4121. The main gas hues are indicated, (b) 
IRS spectrum toward the shell's rim of MB4121. The [Fe II] 24.52 fim line is detected in only one nod. 




sponding column density of CO2, Nco2 — ^ 



Fig. 12.— IRAC 3.6 (top-left), 4.5 (top-right), 5.8 
(bottom-left) and 8.0 /im (bottom-right) images 
(2'x 2') of MB4121. 



wavelength ranges from 13 to 14.7 /im, from 16.2 
to 17.5 /im and from 18.5 to 21 /im in order to 
avoid the emission gas lines of [Co II] 14.74 
and [Fe II] 17.94 fim. A Gaussian is included to 
simulate the blue wing of the 18 fim silicate fea- 



ture. Figure 14(a) illustrate the method and fig- 
ure [T4(b)] shows the resulting optical depth profile 
of the CO2 bending mode. We estimate the corre- 
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10 cm ( Gerakines et al 



To convert this CO2 column density into magni- 
tudes of extinction, one has to k now the abun- 
dance of CO2 relative to hydrogen. IWhittet et al 



(|2007l ) found a tight correlation between Ay and 
Nco2 in the Taurus quiescent intracloud medium, 
Nco2 = qiJ^v - ^0) where q ^ (0.252 ± 0.036) x 
10^^ cm~^ and Aq = 4.3 ±1.0 mag which leads 
to Ay = 27 ± 4 mag toward MB4121. However, 
one shall consider this an upper-limit on the vi- 
sual extinction as the IRDC, which may account 
for most of the CO2 absorption is not represen- 
tative of the matter along the whole line of sight 
toward MB4121. We remove the CO2 absorption 
feature by interpolating the fit of the continuum 
between 14.8 and 16.4 /im. 

Then, our goal is to correct for the silicate 
extinction using the whole SED of the central 
source in MB4121 in a way similar to that used for 
MB3957 (see section l4.3.1.3p . However, in those 
conditions, the fit of MB4121 does not lead to a 
satisfying result: the stellar component fits the 
2MASS and IRAC data but the required extinc- 
tion is too large and the distance too close. We 
thus proceed by step to narrow down the num- 
ber and temperatures of the emission component. 
First, we fit the IRAC and 2MASS observations 
with a single component and extinction. If we 
use a blackbody, then the best fit corresponds to 
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(a) 



(b) 



Fig. 14. — (a) The black solid line is the IRS spectrum toward the central source of MB4121. The red dashed 
line is the polynomial fit to the continuum and the solid blue line is the best fit of the continuum with a 
G aussian included to simula te the blue wing of the 18 /im silicate feature, following the method presented 
in Pontoppidan et al. ( 20081 ). (b) Optical depth of the CO2 bending mode. 



~ 7000 K, IF* - 6 X 10-l^ and rg.y - 4 or 
Ay ~ 30 mag. Assuming the star is a very cool B- 
star, with a radius i?* = 2Rq, its distance is less 
than 20 pc from the Sun, which seems unlikely. 
Moreover, the extinction along the line of sight 
leads to an average density of more than 250 cm~'^, 
about two orders of magnitude higher than the dif- 
fuse ISM. Therefore, we conclude the IRAC and 
2MASS observations reddening is not only due to 
extinction but also to a dusty component. We then 
assume the IRAC and 2MASS observations trace 
an inner envelope of hot dust. We use a v'^B^{T) 
function to fit the IRAC and 2MASS data. The 
best fit is found for Tuot ~ 950 K and Ay ~ 10 
mag. On the other end of the spectrum, we also 
use a i/^Bj^{T) function to fit the IRS spectrum for 
wavelengths A > 19/im, to constrain the proper- 
ties of the continuum originating in the outer shell. 
The best fit corresponds to Tcooi ^ 80 K. Finally, 
we combine the IRAC, 2MASS and IRS observa- 
tions of the central source and fit them with three 
graybody {v'^B„{T)) components: a hot envelope 
of hot dust {That) that accounts for the 2MASS 
and IRAC observations, a cool outer shell (Tcooi) 
that is detected in the MIPS 24 /im images and an 
warm inner component (Ty^arm) that is necessary 
to account for the emission between 10 and 20 /Ltm 
and that could be related to an ejecta. We use our 
preliminary fit of the short and long wavelengths 



SED to put limits on the components tempera- 
tures: 500 < That < 1500 K, 50 < Tcooi < 100 

K and 100 < Ti^arm < 500 K. We run a grid of 
fit for different values of fisM , the contribution of 
the "ISM" extinction curve to the extinction along 
the line of sight, between and 1. The resulting 
depends only slightly on fisM- The best fit 
corresponds to Thot = 1250 K, Tcooi = 74 K and 
Twarm — 296 K and 19 magnitudes of visual ex- 
tinction along the line of sight. The best fit is 
shown in Figure fTSl The agreement with the ob- 
servations is very good, except for the IRAC 3.6 
and 4.5 ^m data points, but these measurements 
might be affected by saturation, as the associated 
flags in the "less reliable" GLIMPSE catalog sug- 
gest. We use the extinction curve derived from 
this fit to correct the MIR gas line fluxes (see Ta- 
ble HU]). 

The visual extinction that we derive from the 
IR SED fitting (19 mag) is in agreement with the 
upper-limit we derive from the CO2 absorption 
feature (27 mag). We obtain a third estimate 
of the extinction thanks to the IRAC 8 /im and 
MIPS 24 /im surface brightness of the IRDC that 
can be seen near MB4121. We extract several 
profiles of the diffuse emission across the IRDC, 
away from point sources, at 8 and 24 /im. The 
surface brightness drops from 80 to 60 MJy/sr 
at 8 /im and from 40 to 30 at 24 /im. Tak- 
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Fig. 15. — Whole IR SED of the central source 
in MB4121. Black diamonds are the 2MASS and 
IRAC data, black solid line is the IRS spectrum. 
The best fit model is in red: dotted and dashed 
lines are the three v^BviT) components, the solid 
line is the total, the diamonds are the transmitted 
fluxes in the 2MASS and IRAC data. 



ing into account the presence of foreground dif- 
fuse emission, we infer the depth at 8 and 24 /im 
due to the IRDC has a lower limit of 0.29 which 
corresponds to Ay > 6 mag, using ^sM 
0.05 and A^a/Av ~ 0.05 (llndebetouw et al 



20051: 



Chapman et al The lower and upper-limits 



of the extinction along the line of sight are signif- 
icantly different but they bracket the value found 
for the best fit of the IR SED. We are therefore 
confident with the correction we use. 

Within this model, we cannot constrain the na- 
ture of the central object since its contribution 
to the NIR SED is not significant. However, we 
put an upper-limit on the dilution factor of 
the star, assuming it does not contribute signif- 
icantly to the 2MASS J band. We assume the 
central source, claimed to be a Be/B[e]/LBV by 
IWachter et al. I (l2010l) . has the radius and temper- 
ature of a B star: R^/Rq = 2 to 7, = 10000 to 
30000 K. We apply the extinction from the best 
fit. We estimate that a B star of 2Rq (resp. 7Rq) 
does not contribute more than 10% to the J band 
if it is farther away than 400 pc (resp. 2.4 kpc). 
However, to remain within the Galaxy, the dis- 
tance to MB4121 cannot be more than a factor of 
a few larger that those estimates. At a distance 
of 400 pc (resp. 2.4 kpc), the MIPS 24 /xm shell 
is about 0.06 pc (resp. 0.36 pc) wide. For com- 



parison, the bright MIPS 24 /xm shell of the LBV 
G79. 26-1-0. 46, at a distance of 1.7 kpc, is about 
0.85 pc wide ( Jimenez-Esteban et al. I2OI0I ). 



4.3.2.2. MIR gas lines The IRS spectrum of 
MB4121 is rich in iron lines, especially toward the 
central source. The spectrum of the outer shell 
shows a weak [Fe II] 25.99 /xm line, and question- 
ably the [Fe II] 24.52 /xm line. The gas line fluxes 
of the lines detected toward the central source 
and corrected for extinction are given in Table [TUl 
These lines are not resolved with IRS {R ^ 600). 
As we show in the previous section, the emission 
toward the central source as seen in the 2MASS, 
IRAC and MIPS 24 /im images, is not dominated 
by the star but by an inner envelope of hot dust. 
We assume the iron lines arise from that envelope 
too. 

We first assume the iron lines fluxes can be ex- 
plained within photoionization models. We use 
the results from the MAPPINGS III model grid 
and compare to the observed brightest iron lines 
([Fe II] 17.94, 24.52 and 25.99 /zm) to constrain 
the temperature of the central source and the den- 
sity of the inner envelope. Figure 16(a) shows 



the [Fe II] 26.0/[Fe II] 24.5 line ratio as a func- 
tion of the [Fe II] 17.9/[Fe II] 26.0 line ratio for 
various gas densities and inner source tempera- 
tures. These iron lines, that we measured with 
a significant signal-to-noise ratio, only lead to a 
rough estimate of the gas density as they all are 
[Fe II] lines. They indicate that, within the range 
of density and temperature that we explore, the 
gas density from where the [Fe II] lines arise has 



to be at least 8 x 10 cm 



The addition of 



dust in the model does not significantly change 
this diagram (see Figure 17(a)). To constrain 
the inner source temperature, we use iron lines 
at different ionization level. In the IRS spec- 
trum of MB4121, the only such line is that of 
[Fe I] 24.04 ^m whose uncertainty is significantly 
larger than that of the [Fe II] lines. Figure [l6(b)| 
shows the [Fe II] 26.0/ [Fe I] 24.0 fine ratio as a 
function of the [Fe 11] 17.9/[Fe II] 26.0 line ra- 
tio for the same range of densities and temper- 
ature as in Figure 16(a) The best match is for 



a 90000 K central source but any temperature 
from 50000 to 110000 K is within l-a of the ob- 
served ratios. The addition of dust to the model 
changes significantly this diagram. For a given 
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[Fe II] 26.0/ [Fe I] 24.0 line ratio, the required 
central source temperature is significantly higher 
when dust is present. Indeed, the photon energies 
at work to ionize the iron atom into Fe+ cations 
are between 7.9 and 16.2 eV (see Table O. The 
dust is competing with the iron at those energy 
level, which corresponds to wavelenghts of 80 to 
160 nm, as dust extin ction is the highest in th e UV 



and far-UV (see e.g. ICompiegne et al.ll201lL and 



reference therein). In MB4121, this competition 
between the extinction by the dust and the ioniza- 
tion of the iron atoms only occurs if the dust is not 
further away from the star than the iron atoms. 
In this case, the observed [Fe II] 26.0/ [Fe I] 24.0 
line ratio is in agreement with any temperature 
above 75000 K. Whether the dust is competing 
with the iron atom or not, the minimum required 
temperatures is about 50000 K which can hardly 
be accounted for by a Be/B[e]/LBV, as Be and 
B[e] stars have an ef fective temperature analog to 
those of B stars fe.g. lLamers et al ]\mi) and LBV 
temperature is wit hin the range 8000 to 40000 K 



( Smith et al. l20o4 ). Within the limits we explore 



with our model, it thus seems difficult to confirm 
the Be/B[e]/LBV na ture of MB4121 suggested by 
Wachter et al. I (l2010l) unless it is an extraordinary 
hot source. 

We then attempt to interpret the observed iron 
line ratios with fast r adiative shocks, using the 
models from I Allen et a l. (2008). A shock is highly 
plausible in the vicinity of a Be/B[e]/LBV star 
as it is likely to undergo mass loss and/or have 
strong winds. Figure [18] shows the same line ra- 
tios as before for models of shock+precursor. We 
use the solar abundance models at two densities 
only, 0.1 and 1000 cm^^, for clarity purposes. In 
both diagrams, the models are shown for a shock 
velocity between 200 and 1000 km.s~^ and for a 
magnetic parameter B/n^^^ between 0.0003 and 
30 fxG.cm^/^. Both diagrams clearly show that 
the low density models cannot account for the ob- 
served ratios, whatever the shock velocity and the 
magnetic param eter. Indeed, only t he highest den- 
sity models from lAllen et al. ( 2008 ) match the ob- 
served line ratios. In both diagrams, the modeled 
curves at 100 cm^"^, not shown here, are distant 
enough from those at 1000 cm~^ that we can rule 
them out. Therefore we conclude that, if the iron 
lines are produced by a fast radiative shock, the 
density in the region from which the iron lines arise 



is 1000 cm~^ within a factor of a few. Without 
better constraint on the density, we do not discuss 
furthermore these diagrams, which do not enable 
us to disentangle the different parameters. Indeed, 
both the shock velocity and magnetic parameter 
required to explain the observed ratios are highly 
dependent on the density. However, they allow 
us to find a set of parameters that would agree 
with the observations: a radiative shock of a few 
100 km.s"^ in a density of ^1000 cm^'^. Since 
none of the lines are resolved with IRS, the veloc- 
ity of a shock in MB4121 might not be higher than 
500 km.s-i. 

We note here that the iron lines arise from 
an envelope near the central source which is un- 
resolved in the MIPS 24 /im image. Therefore, 
its angular size is at most 6". At a distance 
of 400 pc, the minimum distance determined in 
section 14.3.2.11 it corresponds to 2400 AU or 
3.6 X 10^^ cm. With a density of 1000 cm"^^ the 
total column density of the inner envelope is about 
a few 10^^ cm^^ which is equivalent to a visual ex- 
tinction of about 0.02 mag. If MB4121 is farther 
away (e.g. 4 kpc) and the density a factor of a few 
higher than 1000 cm^'', the extinction due to the 
inner envelope is at most 1 mag. 

4.3.2.3. The outer shell The IRS spectra of 
MB4121 both show a continuum at wavelengths 
longer than ~ 20/im. This continuum is detected 
with very similar intensities and spectral shape to- 
ward both position (see Figure 1131) which proves 
it originates in the outer shell of MB4121 detected 
in the MIPS 24 ^m image. In this section, we dis- 
cuss the properties of the dust responsible for the 
continuum emission and compare it with the flux 
measured in the MIPS 70 ^m observations. 

We use the IRS spectrum toward the central 
source of MB4121. The extinction correction per- 
formed on this spectrum (see Figure [15] and Sec- 
tion [13211]) shows that a graybody (v^ x B^{T) 
function) with a temperature of 74 K best fits 
the continuum. We use the dust model from 



Compiegne et all (|2011l) to estimate the dust mass 
in the MIPS 24 /im shell and the mass loss rate 
of the central source. Since we do not have 
any constraint on the dust size distribution in 
the shell, we consider it is dominated by either 
of the three following dust components: small 
amorphous carbon grains (SamC), large amor- 
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Fig. 16.~ Diagrams of (a) the [Fe II] 26.0/[Fe II] 24.5 line ratio and (b) the [Fe II] 26.0/[Fe I] 24.0 fim 
line ratio both as a function of the [Fe II] 17.9/[Fe II] 26.0 fj,m. line ratio, in a dust-free model. The color 
code correspons to (a) the given density (in cm~'^) or (b) the given temperature. All temperatures of the 
inner source are represented by a diamond in (a). Only densities between 6000 and 11000 cm""^, from left to 
right, are represented by a diamond in (b). The cross indicates the observed iron lines ratios after correcting 
for the extinction. 
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Fig. 17. — Same as previous figure but with dusty models. 



phous carbon gr ains (LamC) or amorph ous sil- 
icates (aSil) (see Compiegne et al.l I2OIII for de- 
tails about the component properties). We use a 
20" radius for the MIPS 24 /xm outer sheh. The 
physical radius of the MIPS 24 ^m shell thus is 
about 0.13 X fi*(kpc)pc (i.e. 0.06 and 0.36 pc for 
a 2 and 7 Rq respectively, using the distances 
from section I4.3.2.1|) . The age of the shell then 
is t = 1240 X dH,(kpc)/(wu)mrf/100 km.s~^) years 
(i.e. 500 and 3000 years respectively). To model 
the dust emission we generate the radiation field 
for a B star with a radius between 2 and 7 Rq and 
a temperature set accordingly between 10000 and 



30000 K. We use a distance to MB4121 between 
100 and 5000 pc. We apply to the model spec- 
tra the extinction curve deduced from the whole 
SED fit. We look for the column density of dust 
required to best match the IRS observations for 
wavelength A > 20^m. We find a column density 
of 2 — 3 X 10~^ g.cm~^ if we use SamC only and 
about 1 X 10^^ g.cm"^ if we use LamC or aSil 
only (see Figure [TO)) . We then estimate the total 
mass of dust within the MIPS 24 /xm shell to be 
5.7 X 10^''' X (kpc) M0 if the dust is mainly com- 
posed of large grains and a factor 2 to 3 larger if 
the dust is essentialy made of smaller grains. The 
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Fig. 18.— Diagrams of (a) the [Fe II] 26.0/[Fe II] 24.5 fim line ratio and (b) the [Fe II] 26.0/[Fe I] 24.0 
line ratio both as a function o f the [Fe II] 17.9/[ Fe II] 26.0 /im hne ratio from the MAPPINGS III Library of 



Fast Radiative Shock Models (jAllen et al.ll2008l) . Axis scales are logarithmic. The curves on the left of both 



diagrams are for a density of 0.1 cm while the curves on the right are for a density of 1000 cm ^. The 
red curves are for a given velocity, from 200 to 1000 km.s^^, while the blue curves are for a given magnetic 
parameter B/n^/'^, from 0.0003 to 30 ^G.cm^/^. 



total mass of the shell depends on the dust-to-gas 
mass ratio Xd in the shell that is 1% in the solar 
neighborhood bu t can be as small as 0.25% (e.g. 



Sahai et al. 2000f) . We finally estimate the mass 
loss rate of MB4121 to be 5.7 x lO"'^ x d*(kpc) x 
(0.01/xd) X (w„i„d/100km.s~-^) M0/yr which is in 
agreement with what is expected in LBVs. We fi- 
nally compare the MIPS 70 fj,m flux measured at 
the position of the IRS-LH slit with that extrap- 
olated from the model. The MIPS 70 /im mea- 
sured within the IRS-LH slit toward MB4121 is 
2.0±0.3 Jy, taking into account the uncertainty for 
variations from one "nodding" position to another 
and the photometric uncerta inty of the MIPS - 
GAL 70 ^m observation (15%, ICarev et al.ll2009[) . 
The extrapolated MIPS 70 iim flux from the dust 
model fit with aSil only leads to an extrapolated 
MIPS 70 fim flux of 2.4±0.2 Jy toward the central 
source in agreement with the observed flux. If we 
use SamC or LamC, the extrapolation is more un- 
certain: 3.0±0.4 Jy and 4.1±0.8 Jy, slightly higher 
than the observed value. 

All those estimates are given for the distances 
we determine assuming the central star's contri- 
bution to the 2MASS J band. If the real distance 
is a factor of a few larger than our estimates, the 
timescale to make the shell is lowered by that same 



factor and the dust mass is increased by the square 
of this factor. As a result, the mass loss could be 
higher by a factor of a few. However, the distance 
to MB4121 cannot be significantly higher than a 
few kpc as it surely is a Galactic object. Indeed, 
the modeling of the outer shell continuum emission 
brings an additional constraint on the distance to 
MB4121. As Figure \W\ shows, assuming the tem- 
perature of the central source is not higher than 
30000 K, we can only find a good agreement with 
the observations if the distance to MB4121 is less 
than 8 kpc. If the dust in the MIPS 24 ^m 
shell is dominated by large grains (LamC or aSil) 
then the distance to MB4121 has to be lower than 
^ 2 kpc. The modeling presented here has its 
limitations. For instance, it does only take into 
account radiative heating of the dust due to the 
central source and not collisional heating, which 
would increase the upper-limit on the distance to 
MB4121. 

The combination of 2MASS, IRAC and IRS 
data allows us to model with a limited number of 
components the entire IR SED of MB4121. Con- 
trary to what we model in MB3957, the central 
source of MB4121 is not contributing significantly 
to the IR emission. Therefore, we hardly constrain 
the distance to MB4121 as well as the tempera- 
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ture of its central source. The entire IR SED of 
MB4121 can be explained with three dust com- 
ponents from 75 K in the outer shell to 1250 K 
in the innermost region. Whether there are three 
distinct components or a continuous distribution 
of temperature is beyond the scope of this paper. 
The mod eling of the IRS spectru m with the dust 
model of Compiegne et al.l ( 2011 ) is in agreement 
with the Be/B[e]/LBV nature o f MB4121, as sug- 
gested by IWachter et al. ( 2010[ ). at Galactic dis- 
tances. We note however that the 2MASS, IRAC 
and IRS data were obtained a few years apart from 
each other which may be long enough to observe 



flux variation for LBVs ISmith et al.l ( 2004 ). How- 



ever, within the limits of our modelling, no vari- 
ability is required. 

5. Conclusion 

We present the results from our high-resolution 
Spitzer/IRS mid-IR spectroscopic observations of 
four MBs among the more than 400 detected 
in the MIPSGAL 24 /im survey of the Galactic 
plane. We combine these observations with previ- 
ous 2MASS and IRAC observations, when avail- 
able. We model the mid-IR gas lines detected and 
the dust emission to constrain the physical proper- 
ties of the shell and their central source and there- 
fore determine the nature of these objects. Details 
on each object follow: 

• two MBs (MB4001 and MB4006) exhibit 
a dust-poor IRS spectrum dominated by 
highly-ionized gas lines (e.g. [O IV] 25.9, 
[Ne V] 14.3 and 24.3 fim). Comparison 
with spectra of known planetary nebulae and 
modeling of the gas line ratios with MAP- 
PINGS III leads to the conclusion that these 
two MBs are planetary nebulae with a very 
hot (> 200, OOOK) central white dwarf. We 
relate the almost total absence of dust emis- 
sion features to the evolved stage in which 
the MBs might be. 

• MB3957 exhibits a dust-rich IRS spec- 
trum and has a central source detected in 
the 2MASS and IRAC observations. The 
IRS spectrum is characterized by the pres- 
ence of a few gas lines, among which the 
[Fe III] 22.9 ^m, [S IV] 10.5 ^m and 
[Ne III] 15.6 /im lines which rules out the 



shock excitation mechanism and leads to a 
strong constraint on the inner source tem- 
perature 60000K). Only a white dwarf 
or a Wolf-Rayet star can match such a tem- 
perature. The modeling of the IR SED of 
MB3957, from 1 to 40 /im, enables us to 
constrain the dilution factor of the central 
source and therefore conclude only a Wolf- 
Rayet star is a likely solution. The mod- 
eling also shows that two dust components 
(at 167 and 1750K), about 14 magnitudes 
of visual extinction and a distance of 7.5 
kpc are necessary to interpret the IRS spec- 
trum of MB3957. The MIPS 24 /im dusty 
shell comprises ^ 10~^ M©. Assuming a 
wind velocity of 1000 km/s and a dust-to- 
gas mass ratio of 1%, it implies the outer 
shell of MB3957 is ~ 425 years old and that 
the average mass loss is ~ 10^^ Mo/yr. 

• MB4121 exhibits a dust-rich spectrum and 
has a central source detected in all IR im- 
ages from 2MASS to MIPS 24 /im. Our 
multiple mid-IR spectra toward this object 
reveal many iron lines that arise in the cen- 
tral source and absorption features due to a 
infrared dark cloud along the line of sight. 
It was suggested as a Be/B[e]/LBV thanks 
to near-IR observations of its central source. 
We model the IR emission toward the cen- 
tral source, from 1 to 40 /im with three 
dust components, from 74 to 1250 K and 19 
magnitudes of visual extinction. The cen- 
tral source contribution to the emission is 
required to be small. Therefore, we infer 
MB4121 is at least 400 pc away from us. 
The extinction corrected gas lines cannot be 
fitted by photoionization model but agree 
with fast radiative shock models for density 
of ~ 1000 cm-3. The modelling of the MIPS 
24 /zm outer dust shell results in a total mass 
of a few IQ-^ Mq. It also implies MB4121 
cannot be farther away than 8 kpc. Assum- 
ing a wind velocity of 100 km/s and a dust- 
to-gas mass ratio of 1%, it implies the outer 
sheU of MB4121 is - 1240 years old and that 
the average mass loss is a few 10~^ Mo/yr. 
The distance to MB4121 is uncertain within 
more than an order of magnitude, so are 
these estimates. However, they all fit with 
the Be/B[e]/LBV interpretation. 
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With only four MBs observed in this program, 
we cannot draw any conclusions on the entire set of 
428 objects detected in the Galactic plane. How- 
ever, the combination of ongoing observational 
programs, from the optical to the radio, may soon 
lead to the first statistical results on the MIPS- 
GAL bubbles and the population of evolved mas- 
sive stars in our Galaxy. 
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Table 3: Main parameters used for the MAP- 
PINGS III model 



Main parameters 



Table 5: Gas lines detected in the IRS spectrum of 
MB4006. Units are 10~20 W.cm-^. The lines in 
the SH module may be underestimated by a factor 
ul a few (see text in section 14.2.21 for details) . 



Disable reactions 

Cosmic ray heating 

Geometry 

Radius of Source 

Temperature of Source 

Structure 

Density 

Filling Factor 

Integrate over whole sphere 
Ionization balance 
Photon absorption fraction 
Ionization bounded 



No 
No 

Spherical 
I.OODIO (cm) 
lO'^ - 3 X 10^ K 
Isochoric 

100 - 1.5 X 10^ cm-' 

1.0 

Yes 

Equilibrium 
0.1 

Applies to H 

up to 5% ionization 



Line 


Flux 


Uncertainty 


[S IV] 10.51 ^nJl 


5.05 


0.05 


[Ne II] 12.81 Atm 


0.19 


0.01 


[Ar V] 13.10 Atm 


0.46 


0.02 


[Ne V] 14.32 fj.m 


13.1 


0.9 


[Ne III] 15.56 iim 


3.3 


0.3 


[S III] 18.71 Atm 


1.5 


0.2 


[Ne V] 24.32 ^m 


34.6 


0.3 


[0 IV] 25.89 ^ira 


123.0 


0.5 


[S III] 33.48 fim 


4.9 


0.4 


[Si II] 34.82 fim 


0.83 


0.07 



Dust parameters 



Shattering grain profile 

Power law -3.3 

Graphite grain radius 40-1600 

Silicate grain radius 40-1600 

Graphite grain density 1.8 

Silicate grain density 3.5 
Include PAH molecules 
Fraction of Carbon Dust 

Depletion in PAHs 



Table 6: Line fluxes detected in the IRS spectrum 
ofMB3957. Units are lO'^o W.cm'^. The extinc- 
tion corrected line fluxes results from the whole 
SED fitting (see section 143X31) . 



Line 



Flux Corrected flux 



Yes 

[S IV] 10.51;um 0.9±0.8 4.4±3.9 

0.05 [Ne II] UMum 3.7±0.7 7.0±1.0 

Q- value QHDH < Value PAH photodissociation parameter [Ne III] 15.56/xm 2.3±0.5 5.0±1.0 

PAH switch on Value 1000 [S III] 18.71/im 3.7±0.8 11.0±2.0 

Destroy all C grains with PAHs Yes ^Fe III] 22.92/im 0.4±0.3 1.0±0.9 

[S III] 33.48/im 4.0±1.0 7.0±2.0 

[Si II] 34.82^m 0.9±1.2 1.5±2.0 



Table 4: Gas lines detected in the IRS spectrum of 
MB4001. Units are lO^^o -yv.cm-^. The lines in 
the SH module may be underestimated by a factor 
of a few (see text in section 14.2.21 for details) . 



Line 


Flux 


Uncertainty 


[S IV] 10.51 fim 


2.09 


0.02 


[Ne II] 12.81 urn 


0.015 


0.001 


[Ar V] 13.10 ^m 


0.092 


0.004 


[Ne V] 14.32 ^m 


2.7 


0.2 


[Ne III] 15.56 /im 


0.96 


0.09 


[S III] 18.71 ^ira 


0.49 


0.07 


[Ne V] 24.32 ^m 


6.55 


0.07 


[0 IV] 25.89 fim 


61.7 


0.2 


[S III] 33.48 fim 


1.33 


0.09 



Table 7: Ionization potentials of Iron, Neon, Sulfur 
and Oxygen for their first few ionization states. 
Values are in bold for mid-IR lines that are de- 
tected in MB3957, in italic for mid-IR lines that 
are not detected, and unchanged for ionization 
states that do not emit in the mid-IR. 



Element 


II 


HI 


IV 


V 


Iron 


7.9 


16.19 


30.65 




Sulfur 


10.36 


23.34 


34.79 


47.22 


Oxygen 


13.62 


35.12 


54.93 


77.41 


Neon 


21.56 


40.96 


63.45 


97.12 
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Table 8: IR magnitudes and fluxes of the source at 
the center of MB3957. 



Band 


Magnitude 


Fluxes (mjy) 


J 


13.825±0.034 


4.704±0.156 


H 


11.884±0.039 


18.06±0.665 


K 


10.674±0.023 


35.84±0.792 


3.6 fim 


9.618±0.059 


39.94±2.176 


4.5 fim 


9.041±0.047 


43.48±1.862 


5.8 /xm 


8.732±0.032 


36.96±1.1 


8.0 Aim 


8.490±0.028 


25.78±0.6573 



Table 9: IR magnitudes and fluxes of the source at 
the center of object #4121. 



Band 


Magnitude 


Fluxes (mJy) 


J 


11.711±0.022 


32.97±0.67 


H 


7.889±0.042 


715.7±27.7 


K 


5.82±0.021 


3133±61 


3.6 /xm 


4.915±0.041 


3039±114 


4.5 /xm 


4.519±0.034 


2798±89 


5.8 /zm 


4.046±0.017 


2768±43 


8.0 /xm 


4.054±0.015 


1532±22 



Table 10: Lines fluxes detected in the IRS spec- 
trum of MB4121 toward the central source. Units 
are 10~^° W.cm~^. The extinction corrected line 
fluxes results from the whole SED fitting (see sec- 
tion gXTJl) . 



Table 11: Gas lines detected in the IRS spec- 
trum of MB4121 toward the shell. Units are 
10-20 "W.cm-^. We do not correct those fluxes 
for extinction. 



Line 


Flux 


Uncertainty 


[Fe II] 24.52 fim 


0.1 


0.5 


[Fe II] 25.99 fim 


0.6 


0.3 



Line 


Flux 


Corrected flux 


[Ni II] 10.68/xm 


0.9±0.4 


5.0±2.0 


[He I] or [H I] 12.37Aim 


1.5±0.7 


4.0±2.0 


[Ni II] 12.73Aim 


1.1±0.4 


3.0±1.0 


[Co II] 14.74/xm 


0.3±0.5 


0.8±1.5 


[Fe II] 17.94Aim 


9.0±0.6 


35.0±2.0 


[Fe II] 22.90Aim 


0.3±0.7 


0.9±2.0 


[Fc I] 24.04Atm 


0.3±0.3 


1.0±0.8 


[Fe II] 24.52Aim 


2.6±0.3 


7.3±0.8 


[Fe II] 25.99Aim 


8.4±0.3 


22.0±0.8 


[Si II] 34.82/xm 


4.0±1.0 


7.0±2.0 


[Fe II] 35.35Aim 


3.0±2.0 


4.0±3.0 
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Fig. 19. — contours (solid black lines)for the 
fit of the long wavelength IRS continuum with the 
dust model and infered dust column density con- 
tours (red dashed lines) in /^g.cm~^. The dust size 
distribution is (a) SamC only, (b) LamC only and 
(c) aSil only. 
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